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Abstract. We show that the combinatorial numbers known as Bell numbers are generic in 
quantum physics. This is because they arise in the procedure known as Normal ordering of 
bosons, a procedure which is involved in the evaluation of quantum functions such as the 
canonical partition function of quantum statistical physics, inter alia. In fact, we shall show 
that an evaluation of the non-interacting partition function for a single boson system is identical 
to integrating the exponential generating function of the Bell numbers, which is a device for 
encapsulating a combinatorial sequence in a single function. 

We then introduce a remarkable equality, the Dobinski relation, and use it to indicate why 
renormalisation is necessary in even the simplest of perturbation expansions for a partition 
function. 

Finally we introduce a global algebraic description of this simple model, giving a Hopf 
algebra, which provides a starting point for extensions to more complex physical systems. 



1. Introduction 

The motivation for this work is as follows: Recent analysis of perturbative quantum field theory 
(pQFT) has shown that at the heart of this theory there lies a Hopf algebra description. Since 
pQFT is a subtle and as yet not completely solved area, it is of pedagogical value to examine 
simpler systems to illustrate how these algebraic constructs may arise. In this note we examine 
the simple theory given by the quantum partition function Z of a non-relativistic boson system, 
a completely solvable model (at least in the free boson case). On the basis of a combinatorial 
methodology, we show that certain combinatorial numbers, known as Bell numbers, play a 
central role in such an analysis. The graphical approach to the description of these combinatorial 
numbers leads to a sort of non-relativistic analogue of Feynman diagrams. En route, we also 
remark how that even in this simple case, uncritical expansion methods may lead to divergences, 
which can be cured by regularization; this imitates the renormalization methods of pQFT. This 
theory does not depend on space or time; so it is a non-field theory. However, the description 
of the Quantum Partition Function that we employ uses the coherent state parameter z as 



integration variable; and this plays the role of integration over space and time in the relativistic 
case. 

The core of the note then deals with how an elementary Hopf algebra arises naturally from 
this approach. Although the resulting Hopf algebra, which we call BELL (since it arises directly 
from a graphical description of the combinatorial Bell numbers), is rather too uncomplicated, 
(technically it is commutative and co-commutative), it serves as a starting point for an extension 
to more complicated algebraic systems. Indeed, in later work we show that extensions and 
deformations of this algebra can emulate that of pQFT, thus supplying an alternate, algebraic, 
route from the physical path which takes us from the non-relativistic arena to that of relativistic 
quantum field theory. 



2. Bell and Stirling Numbers; Exponential Generating Function 

The Bell number B(n) simply gives the number of partitions of {1, 2, . . . , n}. Thus for the set 
{1, 2, 3}, we have 

{1,2,3}, {1,2}{3} {1,3}{2} {3,2}{1} {1}{2}{3} 

so that B(3) = 5. Equivalently, B(n) gives the number of ways in which n distinct objects may 
be distributed among n identical containers, some of which may remain empty. The first few 
values of B(n) are 

B(n) = 1, 1, 2, 5, 15, 52, 203, . . . , n = 0, 1, 2, . . . . 

A compact way of defining a combinatorial sequence is by means of a generating function. 
For the Bell numbers, the exponential generating function is given by 
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£5(n)^-=exp(e*-l). (1) 

n=0 

We shall give a proof of Eq.([T]) after we introduce a graphical description for the Bell numbers in 
Section 13.21 A slight generalization of the exponential generating function for the Bell numbers 
is given by that defining the Bell polynomials B n {y): 
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exp (y(e- - 1)) = J2 B n (y) ^ (2) 

With this notation B{n) = B n (l). 

Stirling numbers (of the second kind) S(n, k) are related to the Bell numbers. These are 
defined as the number of ways of putting n different objects into k identical containers, leaving 
none empty. From this definition we have B(n) = J2k=i S( n i 

Bell and Stirling numbers are fundamental in quantum theory. This is because they arise 
naturally in the normal ordering problem of bosonic operators. By the normal ordering operation 
J\ff(a,a^) we mean reorder the boson operators in f(a,a^) so that all annihilation operators 
are on the right. For canonical boson creation and annihilation operators {a^,a} satisfying 
[a, a^] = 1, the Stirling numbers of the second kind S(n, k) intervene through [2J [3] 

n 

M(a^a) n = Y^S{n,k){a)) k a k . (3) 

k=l 

The corresponding Bell numbers B(n) = Ylk=i S( n > ^) are simply the expectation values 

B(n) = {z\(a)a) n \z) z=1 (4) 
taken in the coherent state defined by 

a\z) = z\z) (5) 

for z = 1. 



3. Graphs 

3.1. Bell number diagrams 

We now give a graphical representation of the Bell numbers, based on work of Brody, Bender and 
Meister [HE], which we have extended in [6]. Consider labelled lines which emanate from a white 
dot, the origin, and finish on a black dot, the vertex. We shall allow only one line from each white 
dot but impose no limit on the number of lines ending on a black dot. Clearly this simulates the 
definition of S(n, k) and B(n), with the white dots playing the role of the distinguishable objects, 
whence the lines are labelled, and the black dots that of the indistinguishable containers. The 
identification of the graphs for 1, 2 and 3 lines is given in Figure 1 below. We have concentrated 

n=1: O • 

S(1,1)=1 B(1)=1 



n=2: 




S(2,2)=1 S(2,1)=1 B(2)=2 



n=3: 




S(3,3)=1 S(3,2)=3 S(3,1)=1 B(3)=5 

Figure 1. Graphs for B(n), n = 1, 2, 3. 

on the Bell number sequence and its associated graphs since, as we shall show, there is a sense 
in which this sequence of graphs is generic in the evaluation of the Quantum Partition Function. 
By this we mean the following: The exponential generating function for the Bell polynomials is 
identical to the (integrand of) the Canonical Quantum Partition Function of a non-interacting 
single-mode boson model. Then, when interactions are present these may be incorporated by 
the use of suitable strengths associated with the vertices of the graphs. That is, we can represent 
the perturbation expansion of an interacting model by the same sequence of graphs as in Figure 
1, with suitable vertex multipliers (denoted by the V terms in the same figure). 



3.2. The Connected Graph theorem 

We introduce this important theorem [7] in order to obtain the exponential generating function 
of the Bell numbers. The theorem states: 

Theorem 1 Connected Graph Theorem: If C{x) = Y^n=i c(n)x n /n\ is the exponential 
generating function of a set of CONNECTED labelled graphs then A(x) = expC(x) is the 
exponential generating function of ALL the graphs. 

In Figure [2] we see that there is only one connected graph for each n, so 

C(x) = x/V. + x 2 /2\ + x 3 /3! + . . . = e x - 1 

whence 

A(x) = exp(e x - 1) 
gives the exponential generating function for the Bell numbers. 
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n = 1 n =2 n = 3 
Figure 2. Connected Graphs for B(n), n = 1, 2, 3. 

4. Quantum Partition Function 

We now show that the Canonical Partition Function of Quantum Statistical mechanics is simply 
related to the exponential generating function of the Bell polynomials. 



4-1. Free boson gas and Bell polynomials 

The canonical partition function Z associated with a quantum mechanical hamiltonian H is 
given by 

Z = Trexp(- ( 9#) (6) 

where /3 is the usual inverse temperature ( l/ksoitzmannF)- We start with the elementary case 
of the hamiltonian for the single-mode free boson gas H = ea^a (ignoring an additive constant), 
e > 0. The usual computation of the partition function, exploiting the completeness property 
X)nLo \n)(n\ = I, is immediate: 

Z = Tr exp(— /3ea'a) 
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(l-e-*)\ (7) 



However, we may use any complete set to perform the trace. We choose coherent states as 
defined in Eq.([5]) above, which are explicitly given by 

\ z ) = e ~\ z \ 2 /^Y J {z n /n\)a) n \Q). (8) 

n 

For these states the completeness or resolution of unity property is 

- J d 2 z\z){z\ =1=1 dn(z)\z){z\. (9) 

The appropriate trace calculation is 

Z = - f d 2 z{z\exp(-(3ea j a)\z) = (10) 
vr J 

= - ! d 2 z(z\ : exp(a t a(e- /3e - 1)) : \z) , (11) 
vr J 

where we have used the following well-known relation [HJ [9] for the forgetful normal ordering 
operator '-/(a, aJ) '■ which means normally order the creation and annihilation operators in f 
forgetting^ the commutation relation [a, a)] = 1: 

Nexp{xa)a) =: exp(a t a(e x - 1)) : . (12) 

1 This procedure generally alters the value of the operator to which it is applied. 



We therefore obtain, integrating over the angle variable 9 and the radial variable r = \z\, 
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Z=-J d6 J rdrexp (r 2 (e~ pe - 1)) , (13) 

which gives us Z = (1 — e - ^) -1 as before. 

We rewrite the above equation to show the connection with our previously-defined 
combinatorial numbers. Writing y = r 2 and x = —fie., Eq. (jl3p becomes 

oo 

Z = J dyexv(y(e x -l)). (14) 
o 

This is an integral over the classical exponential generating function for the Bell polynomials as 
given in Eq.([2]). This leads to the combinatorial form for the partition function 
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Z= / dyY^B n (y)^-. (15) 

4-2. General partition functions 

We now apply this graphical approach to the general partition function in second quantized 
form, with the usual definition for the partition function Eq.Q. In general the hamiltonian is 
given by H = ew(a, a^), where e is the energy scale, and w is a string (= sum of products of 
positive powers) of boson creation and annihilation operators. The partition function integrand 
F for which we seek to give a graphical expansion, is 



Z(x) = J F(x,z)dfi{z) , (16) 

where 

F(x,z) = (z\ exp(xw)\z) = (x = — j3e) 
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(17) 



with obvious definitions of W n and V n . The sequences {W n } and {V^} may each be recursively 
obtained from the other [10 . This relates the sequence of multipliers {V^} of Figure 1 to 
the hamiltonian of Eq.(j6j). The lower limit 1 in the V n summation is a consequence of the 
normalization of the coherent state \z). 

Although Eq. (|15p is remarkably simple in form, it is often by no means a straightforward 
matter to evaluate the analogous integral for other than the free boson system considered here. 
Further, as we show below, we may not interchange the integration and the summation, as each 
individual y integral diverges. At this point the analogue of this model to perturbative quantum 
field theory becomes more transparent. For example, to calculate the partition Function Z(x), 
we must integrate over the Partition Function Integrand F(x,z). Here, the integration over 
the coherent state parameter z plays the role of the space-time integration of Quantum Field 



Theory. As we have just remarked, with reference to Eg. (|15j) for the non-interacting case, and 
as in pQFT (perturbative quantum field theory), term-by-term integration, corresponding to 
fallaciously interchanging the summation and integration actions, results in infinities at each 
term, as we now verify for the non-interacting case. Thus even here an elementary form of 
renormalization is necessary, if we resort to term-wise integration. 

Considerations such as these have led many authors, as noted in reference [I], to consider a 
global, algebraic approach to pQFT; and this we shall do in Section [UJ 

But first we show how divergences arise. 



5. The Dobiriski Relation and Divergence 

5.1. Divergence 

The celebrated Dobihski relation for the Bell numbers is given by (see, for example, 

1 00 k n 

5 ( n ) = -E7T' " = 1,2,... (18) 
e k=o K 



and for the Bell polynomials by 
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B n {y) = e"^-/. (19) 
fc=o K - 

These relations have been extended in [12j . 

Using the Dobihski relation Eq. (|19p we can rewrite the expression for the partition function 
for the free single- mode boson Eq. (|15p as 

°o n 

dyY,Bn(y)^ (20) 
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We may attempt to write this as a perturbation expansion (in x = —fie) by interchanging the 
integration and summation, thus: 

n=0 k=0 

n=0 k=0 

oo n oo 

= E^rE**- ( 24 ) 

n=0 fc=0 

Therefore the resulting "perturbation expansion" is of the form 
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n - k=0 



with each term in the expansion being infinite. 

Thus even in this simple case of a free boson expansion, a perturbation expansion obtained 
from the illegal interchange of integration and summation is divergent in every term. 



5.2. Regularization 

In this section we illustrate how the method of regularization may be applied to alleviate the 
problem of divergences. Writing the partition function Eq. (|14p as 



Z = / e~ ay dy (a = 1 - e _/Je ) (25) 



immediate evaluation gives Z = 1/a as in Eq.©. 
However, if we expand 

Z= J2 (-<*y) n /n\ dy (26) 

J ° n=0 

and exchange the order of integration, we obtain 



OO /"OO 

(FALSE) Z=J2 {-cc) n /n\ / y n dy (27) 
ZTn Jo 



n=0 

each term of which diverges. The regularization method is to write 

i-M 00 

Z =^ m / E(-«y) n /^^ (28) 

and, since the integration is over a finite range, we may interchange integration and summation, 
to obtain 

»M 
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Z = lim Y {-a) n /n\ / y n dy (29) 
M ->°°^o Jo 

oo 

lim V{(-a)7n!}{M( n+1 V(n + l)} (30) 



M-5>oo 

n=0 

oo 

lim y{(-aM)^/(n + l)!}(-l/a) (31) 

n=0 



-Ma 

1/a (33) 



lim {(e- Ma -l)(-l/a)} (32) 

Ai— >oo 



as before. 

Note that y bounded by a finite M implies the eminently reasonable assumption of a maximum 
energy for the coherent state as defined in Eq.©. 



6. Hopf Algebra structure 

So far we have discussed the analytic and combinatorial structure of our simple model. Modern 
approaches to perturbative quantum field theory introduce a global algebraic structure. We 
describe the analogue of this Hopf Algebra structure in the context of our model, which we refer 
to as BELL below. We first introduce a basic Hopf Algebra, which we call POLY, generated by 
a single parameter x. This is essentially a standard one- variable polynomial algebra on which we 
impose the Hopf operation of co-product, co-unit and antipode, as a useful pedagogical device 
for describing their properties. The Hopf structure associated with our model described above 
is simply a multi-variable extension of POLY. 



6.1. The Hopf Algebra POLY 

POLY consists of polynomials in x (say, over the real field 1Z, for example). The standard algebra 
structure of addition and associative multiplication is obtained in the usual way, by polynomial 
addition and multiplication. The additional Hopf operations are: 

(i) The coproduct A : POLY — ^POLY x POLY is defined by 

A(e) = e x e (unit e) 
A(x) = x x e + e x x (generator x) 
A(AB) = A(A)A(B) otherwise 

so that A is an algebra homomorphism. 

(ii) The co-unit e satisfies e(e) = 1 otherwise e(A) = 0. 

(iii) The antipode S : POLY — ^POLY satisfies 5(e) = e; on the generator x, S(x) = —x. It 
is an anti-homomorphism, i.e. S(AB) = S(B)S(A). 

It may be shown that the foregoing structure POLY satisfies the axioms of a commutative, 
co-commutative Hopf algebra. 

6.2. BELL 

We now briefly describe the Hopf algebra BELL which is appropriate for the diagram structure 
introduced in this note, defined by the diagrams of Figure 1. 

(i) Each distinct diagram is an individual basis element of BELL; thus the dimension is infinite. 
(Visualise each diagram in a "box".) The sum of two diagrams is simply the two boxes 
containing the diagrams. Scalar multiples are formal; for example, they may be provided 
by the V coefficients. 

(ii) The identity element e is the empty diagram (an empty box). 

(iii) Multiplication is the juxtaposition of two diagrams within the same "box". BELL is 
generated by the connected diagrams; this is a consequence of the Connected Graph 
Theorem (as discussed above in 13. 2p . Since we have not here specified an order for the 
juxtaposition, multiplication is commutative. 

(iv) The coproduct A : BELL — >BELL x BELL is defined by 

A(e) = e x e (unit e) 
A(y) = yxe + exy (generator y) 
A(AB) = A(A)A(B) otherwise 

so that A is an algebra homomorphism. 

(v) The co-unit e satisfies e(e) = 1 otherwise e(A) = 0. 

(vi) The antipode S : BELL — s>BELL satisfies S(e) = e; on a generator x, S(y) = —y. It is 
an anti-homomorphism, i.e. S(AB) = S(B)S(A). 

6.3. BELL as an extension q/POLY 

It may be seen that BELL is a multivariable version of POLY. To show this, we code the 
diagrams by letters. We use an infinite alphabet Y = {2/1,2/2,2/3 • • •} = {yk}k>i and code each 
connected diagram with one black spot and k white spots with the letter y^. An unconnected 
diagram will be coded by the product of its letters. 
In this way, the diagrams of Figure 1 are coded as follows: 



• first line : y\ 

• second line : y\ and yi 

• third line : y\ and y\yi and 7/3 . 

Thus one sees that each diagram of weight n with k connected components is coded bijectively 
by a monomial of weight n (the weight of a monomial y% x y%^ ■ ■ ■ y% r is just the sum of the indices 
Y^j=\ ij) an< i k letters. The algebra BELL is coded by commutative polynomials in the infinite 
alphabet Y; that is, the coding is an isomorphism BELL — > 7Z[Y]. As an aside, one may note 
that the basis elements of BELL are sometimes referred to as forests. As above, it may be shown 
that the foregoing structure BELL satisfies the axioms of a commutative, co-commutative Hopf 
algebra. 

7. Discussion 

In this note we have shown the value of the combinatorial method, which leads inter alia through 
the graphical approach, to a Hopf algebra description of a simple non-relativistic free boson 
system. We named the resulting (commutative, co-commutative) Hopf algebra BELL. To extend 
this simple description to more realistic models, it is convenient to introduce a symmetrized 
form of BELL; i.e. where the black and white spots appear symmetrically [13]. Thus we have 
introduced the Hopf algebra DIAG where a typical series of graphs is given in Figure [3] below. 
Note that since the black and white spots play symmetrical roles, in the figure we permit 
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Figure 3. Some Graphs of DIAG 

ourselves to exchange their order from that used in, for example, Figure [TJ We may also apply 
"weights" L m to the black spots as well. Although the resulting Hopf algebra DIAG is still 
commutative/co-commutative, it may in turn be extended to a non-commutative Hopf algebra 
(LDIAG) which forms the basis of an extension of this algebraic approach to a description 
sufficiently general to emulate those algebras which describe pQFT. This extension is described 
in detail in a subsequent artcle in this conference [14] . 
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